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Abstract. A finite array of N globally coupled Stratonovich models exhibits a 
continuous nonequilibrium phase transition. In the limit of strong coupling there 
is a clear separation of time scales of center of mass and relative coordinates. The 
latter relax very fast to zero and the array behaves as a single entity described by 
the center of mass coordinate. We compute analytically the stationary probability 
and the moments of the center of mass coordinate. The scaling behaviour of the 
moments near the critical value of the control parameter ac{N) is determined. 
We identify a crossover from linear to square root scaling with increasing distance 
from Og. The crossover point approaches Oc in the limit N oo which reproduces 
previous results for infinite arrays. The results arc obtained in both the Fokker- 
Planck and the Langevin approach and are corroborated by numerical simulations. 
For a general class of models we show that the transition manifold in the parameter 
space depends on N and is determined by the scaling behaviour near a fixed point 
of the stochastic flow. 
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1. Introduction 



Arrays of stochastically driven nonlinear dynamical systems may exhibit 
nonequihbrium phase transitions of continuous or discontinuous type, for a recent 
review see P, cf. also [21 [3]. Concepts developed to describe equihbrium phase 
transitions such as symmetry or ergodicity breaking, order parameter, critical 
behaviour, critical exponents etc. have been successfully transfered to noise induced 
nonequihbrium phase transitions. The structure of the theory will be generically 
of mean field type, if infinite globally coupled arrays are studied which allows for a 
number of analytical results. 

Remarkably, essential characteristics of phase transitions can already be found 
in the case of a single Stratonovich model. This is mainly due to the multiplicative 
nature of the noise. Models driven by additive noise do not show this pecuhar 
property. The Langevin equation for the single-site Stratonovich model [H El E] 
reads 

dx = {ax - x^)dt + ax o dW{t), (1) 

where a is a control parameter, a denotes the strength of the noise and W{t) is 
a Wiener process with autocorrelation {W{t)W{s)) = min(t,s). Equation Q is 
interpreted in the Stratonovich sense as indicated by the symbol o. The Stratonovich 
model describes, e.g., the overdamped motion in a biquadratic potential U{x) = 
— |a;^ + where the control parameter is stochastically modulated, a — > a + 
with a Gaussian white noise 

The associated Fokker-Planck equation (FPE) describing the evolution of the 
probability density P{x,t) is 

dtP = -d,{[{a - ^)x -x'- yX^Sjp}. (2) 

Equation ([2| has a weak stationary solution, a Dirac distribution 5{x) located at 
the common zero x = of drift and diffusion coefficient, which is also a zero of the 
stochastic fiow in Eq. ([T]). If the system is initially at x = it wiU always stay there. 

Furthermore, there exist spatially extended strong stationary solutions 
determined up to a constant factor, Ps{x) oc exp{ — (x/a)^}. Ps{x) wiU live 

on 5"+ = [0, oo) if the initial distribution lives on S'+\0, and on S- = (— oo,0], if the 
initial distribution lives on S'_\0. The constant is determined such that the solution 
is normalized integrating over the support and can be interpreted as probability 
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density, i.e. 

Ps{x) = \\xr'"'-'e-^^'"'^\ (3) 
Z = I rfx|x|2'^/'^'-^e-("/'^)' = -a2"/<^'r(a/(T2), (4) 

provided 2a/ a'^ > 0. For 2a/ a'^ < the normalization Z diverges since the integrand 
in Q scales hke as a; — > 0. In this case it can be shown (see Appendix 



[Bj) that a weakly normalized version converges to the known weak solution 6{x) and 
X = is an absorbing fixed point of the system. 

If fractions of the initial distribution of given weights live on S-, on 5*+, and on 
0, all will keep their weight and evolve to the stationary probability densities living 
on their respective support as guaranteed by a if-theorem [7]. 

The Stratonovich model exhibits a strong ergodicity breaking |8] depending on 
the control parameter a, since the state space decomposes into regions where the 
system cannot reach one region if it has started in a different one. For a < the only 
stationary solution is 5{x), i.e. the fixed point x = of the stochastic dynamics is 
absorbing. Additionally, for a > we have the spatially extended solution (|3j living 
on S± depending on the initial distribution. This is refiected by the mean value 

fo ifa<0, 

(x)± = / dxxPsix) = \ T{a/a^ + 1/2) ^ (5) 

Js± \ ±a ^ ' ' ifa>0. 

Obviously, {x)± can serve as an order parameter and shows a critical behavior 

(x)± ~ ±^(a - ac{l)f as a ^ 0^(1) = with (3 = 1. 

Note that also the location of the maximum of the spatially extended density 

undergoes a bifurcation, x™^^ = for < a < a™^^(l) = a'^/2 and x™^^ = 
±(a - a™^)i/2 for ^ > amax(^^)_ 

The critical behaviour of an array of infinitely many globally coupled 
Stratonovich models has been thoroughly investigated in |9]. The scahng of 
higher moments was considered in [10], see also [H]. The stationary probability 
density is the solution of a nonlinear Fokker-Planck equation which depends on 
the order parameter. The scaling behaviour of the order parameter is analytically 
determined, {x)± ~ ±(a — ac{oo)Y as a — ac(oo) with ac(oo) = —a^/2 and 
(3 = sup{l/2, cT^/(2D)}, where D is the strength of the harmonic couphng between 
the systems |i9j. The strong coupling limit, D oo, of an infinite array of globally 
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coupled systems was analytically treated already in the pioneering paper [T2|, cf. 
also pllll]. 

In this paper we investigate nonequilibrium phase transitions in finite arrays 
of globally coupled Stratonovich models in the strong coupling limit. We introduce 
center of mass and relative coordinates and exploit that for strong couphng there is 
a clear separation of time scales. The relative coordinates relax very quickly to zero 
and the system behaves as a single entity described by the center of mass coordinate 
Rt. Thus, we can adiabatically eliminate the relative coordinates. The stationary 
probability density of the center of mass coordinate Ps{R) is analytically calculated 
for a class of nonlinear systems and a scheme to determine the transition manifold in 
the parameter space is developed. For finite arrays of Stratonovich models the mean 
value {R) of the center of mass coordinate is computed analytically. Near a critical 
value of the control parameter a the stochastic system shows a scaling behaviour 
similar to the order parameter of the single Stratonovich model with the same critical 
exponent /? = 1 but with a different ac{N) which is also given analytically. Keeping 
a finite small distance to adN) we recover for — > cxd the known result of the self 
consistent theory [9] with critical exponent (3 = 1/2, see above. For finite we 
identify a crossover value of the control parameter a^{N). For adN) < a <^ o-i^iN) 
we have a linear scaling as for = 1 whereas for a ^ a^(A^) a square root behaviour 
as for ^ oo is observed. The analytical results are coroborrated by numerical 
simulations. 

Recently, finite arrays of (non-) linear stochastic systems have been investigated 
also by Mufioz et al. ^0\, and by Hasegawa [T5] . 

Mufioz et al. tried to obtain for multiplicative noise characteristics of the 
probability density of the mean field for finite A^. They argued that the Langevin 
equation for the mean field variable is of similar form as the Langevin equation for 
a single system. Assuming that the multiplicative driving noise and the local field 
variable are uncorrelated, they inferred the scaling behaviour of the variance of an 
effective multiplicative noise with A^, and of the critical value adN) of the control 
parameter. They also predicted a crossover from a critical exponent /? = 1 near 
ac{N) to the critical exponents for A^ — > oo for larger distances to adN). Note that 
in lin] the Langevin equation was treated in the Ito-sense which leads to a shift of the 
critical control parameter compared to the same equation in the Stratonovich-sense. 

Hasegawa considered finite systems with additive and multiplicative noise using 
his augmented moment method which is applicable for small noise strength. He 
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emphasized that multiphcative noise and the local field variable are not uncorrelated 
in contrast to the assumption in [lO] and demonstrated some consequences of such 
a simplification. 

Our approach, though similar in spirit to [lO], is controllable, valid in leading 
order for strong couphng D, and provides exphcit analytic results which are confirmed 
by independent numerical simulations. It may serve as a starting point to calculate 
next order corrections ~ 

The paper is organized as follows. In the next section we consider two 
harmonically coupled Stratonovich models and show that for strong couphng the 
center of mass coordinate R is the relevant degree of freedom. The mean value 
of R shows a critical behavior which is analytically characterized. Section [3] deals 
with a class of globally coupled systems of general kind. For strong coupling we 
compute analytically the stationary probability distribution Ps(-R) after eliminating 
the relative coordinates. Further, we determine the transition manifold in the 
parameter space where Ps{R) undergoes a transition from a delta-distribution to a 
spatially extended solution. In Sec. [4]we specialize to the case of globally coupled 
Stratonovich models and determine the critical behaviour of the order parameter 
and of higher moments of R for strong coupling. Conclusions are drawn and a 
summary is given in Section [5| In [Appendix A| we introduce the concept of weak 



normalization for the case that a spatially extended solution of the stationary FPE 



cannot be normalized in the naive sense. Appendix B shows that the Langevin 
approach both in Stratonovich- and in Ito-interpretation leads to the same results as 
the Fokker-Planck approach used in the main part of the paper. 

2, Two coupled Stratonovich systems 

We consider a pair of particles with coordinates xi{t) and X2(t) in a biquadratic 
potential which are coupled harmonically with positive coupling strength D and 
each subjected to independent Gaussian white noise of strength a. The system of 
Langevin equations reads 

„3 



axi — Xi — D {xi — Xj) dt + axi o dWi(t), i = 1,2 , (6) 

where Wi{t) denotes independent Wiener processes with {Wi{t)Wj{s)) = 
6ij min {t, s). In contrast to Eq. ([T]) no exact solution of system (|6]) is known. 
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The joint probability density P{xi,X2;t) is governed by the FPE 



=1,2 



=1,2 



(7) 



where, adopting the notation of [21] . 

Di{xi, X2) = (a + —)xi -x] - D - x^), 



A,i(a;i) 



cr 



2 -^j^u 



(8) 
(9) 



denote drift and diffusion coefficients, respectively. 

One can show that the system ([T]) exhibits no detailed balance. Hence, there is 
no easy way to obtain analytically the stationary solution P,- 

For strong couphng, however, a systematic analytical approach is possible. With 
increasing coupling strength the particles become tightly glued together and move as 
a single entity. Therefore it appears natural to introduce center of mass and relative 
coordinates. Simulations of Eq. ([6| show that indeed the stationary distribution of 
the relative coordinate ps{r) becomes very sharp for large values of D, cf. Fig. [!} 




Figure 1: Stationary probability densities for center of mass Ps{R) (left) and relative coordinates 
Ps{r) (right) for two coupled systems. The distribution of relative coordinates is symmetric with 
respect to zero and becomes very sharp with increasing strength of the coupling D. The symbols 
show histograms from 4 x 10^ realizations obtained by solving Eq. ^ with a stochastic Runge- 
Kutta scheme |16]- Parameters are D = 1 (squares), 10 (circles), and 100 (triangles); a = 1 and 
cr^ = 1. Initial values were all in the positive sector. Entries of several bins are omitted to avoid 
overloading; the lines are guides to the eye. 
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The stationary distribution of the center of mass Ps{R) shows a behaviour which 
is similar to the distribution of a single Stratonovich model. For large values of a 
we have a monomodal distribution which vanishes at the boundaries of the support, 
cf. Fig. [T| For small values of a the distribution Ps{R) diverges as i? — in 
a normalizable way, cf. Fig. [2} For even smaller values of a all trajectories Xi{t) 
approach zero and the distributions of both r and R are (^-distributions. Accordingly, 
the mean value (R) undergoes a continuous transition at a critical value of a. In the 
following we analytically calculate Ps(-R) and {R) and its scaling characteristics in 
the strong couphng limit, D — > oo. We introduce the center of mass coordinate R{t) 




Figure 2: Stationary probability density for center of mass Ps{R) (left) and relative coordinates 
Ps{r) (right) for two coupled systems in a semilogaritlimic plot. For sufficiently small control 
parameter, Ps{R) diverges for i? ^ in a normalizable way. Ps{r) becomes essentially sharper 
compared to Fig. [l] Parameters are a = —0.05, D = 100, ~ I (triangles) and 0.49 (circles). The 
symbols show data from 5 x 10^ realizations generated by a stochastic Runge-Kutta algorithm |16] . 
Initial values were chosen in the positive sector. The lines are guides to the eye. 

and the relative coordinate r(t) by 

R=^{xi+X2), r = ^{xi-X2) (10) 
with the inverse transformation 

Xi = R + r, X2 = R — r. (11) 



With 



d.,^, = ^idn±dr), (12) 

dl,, = \{dl + d'^±dl), (13) 



the Langevin equations ^ then transform to 
dR 



a 



{aR -R^- 3Rr^) dt + — [R o dWi{t) + r o dW2{t) 

v2 



[U) 



dr = [{a-2D)r-r^-3rR^]dt 



a 

71 



where the transformed Wiener processes Wi{t) are defined as 



Wi 
W2 



1 

71 



1 1 
1 -1 



Wi 

W2 



r o dWi{t) + R o dW2{t)] , (15) 



(16) 



with {Wi{t)Wj{s)) = 5ijYmn{t,s). 

The FPE associated to ( ppSj ) governing the probabihty density of center of 
mass and relative coordinates P(-R, r; t) reads 



dtP = {Cr + Cr + Cru) P.. 

where the Fokker-Planck operators are 
= - Or 



2 2 
+ y )/? -R"- 3Rr^ - ^Or {R^ + r^) 



dr 



[a — 2D + —)r — r" 



3rR^-^dr (R^ + r^) 



o^dl^ rR . 



(17) 

(18) 

(19) 
(20) 



Note that only Lr depends on D. In the strong couphng hmit D ^ 00 the relative 
coordinate vanishes, — > 0, on a very fast time scale of the order 1/D, cf. Eq. (15). 
Hence, the stationary probability density factorizes to Ps{R,r) = ps{R)6{r) with a 
Dirac distribution for the relative coordinate. In this case there is no flow related to 
the relative coordinate r, i.e. Cr-P = (^trP = 0, since for any suitable function ip 



drdr\'p{r)5{r)\ = 



Integrating Eq. with respect to r yields in the stationary case 

: 



drCRPs = -d^ 



'R 



(aR - i?3 - ^R^dR)ps 



(21) 



(22) 



Similarly as for the single Stratonovich model, there is always a weak solution 5{R). 

For initial values Xi{Q) > Vi (or Xi{0) < Vi) the spatially extended solution 
of Eq. (22) lives on the support 5*+ = [0, 00) (or on = (— oo,0]) and can be 



normalized provided a > ac(2) = — o"^/4. For a < ac{2) the weakly normalized 
version of the spatially extended solution converges to 6{R). Thus, we have 

5{R) for a < 0^(2), 

Ps{R) = < 1, „ (23) 



|^|4a/.2g-2fiV-^ for a > a,(2), 



^2a/.2+i/2p^2a/a2 + l/2). (24) 
There is a strong ergodicity breaking when a crosses ac(2). The mean value {R)± 



calculated with (23) is 



{Ry_ 



ifa<ac(2), 

ifa>ac(2). 



(25) 



"v^r(2a/a2 + l/2) 
and scales like {R)± ~ ±^(a - ac(2))^ with ^ = 1 as a 0^(2) . 

3. General N-site systems 

3.1. Adiabatic elimination of relative coordinates 

In the following we demonstrate that the strategy sketched above can be generalized 
for a class of N coupled systems. We consider 

N 



D 



N 



Y^^ixi-Xj) dt + g{xi)odWi{t), 



(26) 



with i = 1, . . . ,N and where / and g are smooth (with no singularities) and twice 
differentiable chosen such that the stochastic process x(t) = {xi{t),i = 1,...,A^} 
has natural boundaries at infinity |H [20]. Both / and g may depend on a d- 
dimensional set of control parameters a. D > is the couphng strength of the 
harmonic attraction. Note that we have absorbed a factor a, the strength of the 
noise, in the function g. 

The FPE for the joint probability density P(x;t), x = {xi,i = 1,...,A^}, 



associated to (26) reads 



N N 

1=1 j=i 



(27) 
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Using the shorthands fi = f{xi), gi = g{xi), and = d^^gi, drift coefficient and 
diffusion matrix are given by 

1 D ^ 

A = fi + ^g'iQi - Yl " ' (2^) 

D^,J = l9^S^,. (29) 

It is advantageous to introduce center of mass and relative coordinates {R, r}, 
r = {rfc, k = 2, . . . , N}, by the hnear transformations 

1 ^ 

^ = ]^E^- (30) 
j=i 

rk = Xk- R for = 2, . . . , iV . (31) 

The inverse transformation is given by 

TV 

xi=R-Y,n, (32) 

k=2 

Xk = R + rk for k = 2,...,N . (33) 



Observing the rules for linear transformations we have 

+ (34) 

i=l * A;=2 
i,j=l ■' k=2 k,l=2 

Drift and diffusion coefficients in the new coordinates are given by, cf. also [2T| . 

«.-EgA. I'.^Eg^A.. (36) 

i=l *ii=l 

where and 2; stand for the new coordinates R, r^, and r^, respectively. 

Again, the FPE determining P{R,r;t) has the form dtP = CP, with C = 
Cr + Cr + CrR where 

CR=-dR{DR-dRDR^R), (37) 

N N 

l^r = -J29rdDr,-J29r,Dr„r,), (38) 
k=2 1=2 
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N 

^rR = J2<.^^^r,. (39) 

k=2 

Explicitly, the new drift and diffusion coefficients are 

1 ^ 1 

^« = {f^ + 2^^9^) ^ (40) 

i=l 

1 

Dr, = -Dr + fk + -g'^Qk - D——rk , (41) 



2"^" iV 

TV 



2iV2 

i=l 



^iJ,r, = Afe.R = ^^fc - Dr^r , (43) 

^r-„n = Dr^r -^{gl + gf) + l^^^ _ (44) 

All arguments in /j, ^(j, and have to be expressed by {R, r), see Eqs. ( 32p3 ). Note 
that only Dr^^ depends on the coupling strength D explicitly. 

For large times the probability density P{R, r; t) converges to a stationary 
probability density, cf. [7], determined by CPs{R,r) = 0. For D oo this enforces 

5^— [r,P,(i?,r)] = 0, (45) 

k=2 ^ 

which has a weak solution 

Ps{R,r)=ps{R)6{r). (46) 

In the strong couphng limit all fluctuations of the relative coordinates vanish. The 
system is concentrated on the center of mass and moves stochastically as a whole, 
combined particle. 

The probability density of the center of mass Ps{R) can be determined by 
integrating Ps{R,r) over all relative coordinates. Performing this integration we 
obtain from the stationary FPE 

j d^'-h CPs = j rf^"^r Cr{R, r)PsiR, r) = 0, (47) 

provided that the boundary terms associated with the relative coordinates vanish. 
In the strong coupling hmit we have Ps{R, r) oc 5(r), (47) holds in any case and leads 
to 

CRiR,O)psiR)=0 , (48) 
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where Cr is given by (37). From (40) and (42) we infer drift and diffusion for r = 

as 

DR{R,Q) = f{R) + ]^g\R)g{R), (49) 
DR,R{R,0) = ^g\R). (50) 



The spatially extended strong solution of (48) is given by 



Ps{R) = \\9{Rr-' ^y^v{^Nj dR!^^] (51) 

provided that the normalization constant Z is finite. Whether or not this is the case 
depends on the scahng behaviour of the functions f{R) and g{R) near a common 
zero i?o which, if existing, will build a boundary of the support. This is explained in 
detail in the next subsection. 



Equation (50) shows that in the strong coupling limit the diffusion coefficient 
Dufi scales like cr^/A^, cf. also jTOl [I5]. For the infinite system and finite noise 
strength a the stationary probability density of the center of mass Ps{R) is a Dirac 



measure located at one of the attractive zeros of the drift coefficient (49), depending 
on the initial conditions. 

For D ^ oo all particles are strongly correlated. The variance of the coordinate 
Xi{R,r) of an arbitrary system i calculated with Ps{R,r) = ps{R)5{r) is 

(x?) - {x,f = {R') - {R)'. (52) 

Due to the strong correlations, the variance of the center of mass scales hke iV° in 
contrast to the case of uncorrelated systems where the central limit theorem predicts 
a scahng hke N~^. 



3.2. Determination of the transition manifold 

There will be a strong ergodicity breaking if the state space decomposes into different 
regions with the property that one region will not be accessible if we start in a 
different one |B]. 

For multiplicative noise, zeros of the stochastic flow separate the state space into 
mutually non-accessible regions. If we place the system initially on such a zero, i.e. 
on a fixed point of the stochastic dynamics, it will stay there forever. Accordingly, 
the FPE has a weak solution, a 5-distribution living on that fixed point. If any 
trajectory in the neighborhood 'asymptotically' reaches the fixed point (the fixed 
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point is absorbing), there will be no spatially extended probability density in this 
neighborhood. The spatially extended stationary solution of the FPE cannot be 
normalized in the naive sense. The weak normalization procedure leads to the weak 
solution. 

If trajectories cannot reach the fixed point, the stationary solution of the FPE 
will be normalizable and we will have a spatially extended probability density living 
on the support bounded by the fixed point. This properties can be exploited to 
determine the transition manifold in the parameter space. 

We suppose that /(i?, a) and g{R, a) near a common zero Rq have the following 
scaling behaviour 

f{Ro + e)r^Afe^f, (53) 
g{Ro + e) ^Ags^^, (54) 



where mf,mg > 0. Near _Ro we have for the stationary solution (51) of the reduced 
FPE (|4§ 

Ro+e 

PsiRo + s) oc \er^(^-^)exp{2Nj dR'j^^}. (55) 

For nif — 2mg > — 1 the integral in (55) gives a contribution oc 5™/-2mg+i 
the upper boundary which vanishes for e — > so that in leading order Ps{Ro + £) oc 
\£\™-3('^-'^\ The exponent mg{N — 2) is negative only for = 1. In this case, if 
Trig ^ 1 the singularity of ps at R = Rq is not normalizable and we have only a weak 
stationary solution Ps{R) = 6{R — Rq). Note that for ^ 2 coupled systems of this 
kind the singularity of Ps does not occur. 

For TTif — 2mg = —1 the integral gives a logarithmic contribution (Af/A^) In \e\ 
which leads to 

Ps{Ro + e) oc \e\"^AN-2n2NAf/Al (gg) 

If the exponent in (56) is smaller than —1 the density Ps{R) will diverge for R ^ Rq 
and will not be normalizable in a naive way. The weak normalization procedure 
leads to a Dirac measure located at Rq. If the exponent is larger than —1 the density 
Ps{R) will be normalizable and we will have a spatially extended probability density. 
The transition manifold Ac in the control parameter space M*^ is determined by the 
condition that the exponent in (56) is —1, 

A = {a G M'^ : mg{N - 2) + 2NAf/Al = -l} . (57) 
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For nif — 2mg < —1 the integral in (55) diverges like —Af / {A^lrrif — 2mg + 
as e ^ 0. Accordingly, i?o) = for A/ > and 

Ps{R — > Ro) = oo for Af < 0. In the first case Ps(-R) is normalizable and we 
have a spatially extended stationary probability density. In the latter case the weak 
normalization procedure yields a Dirac measure at i?o- A change in the sign of Af 
induced by tuning a control parameter is associated with a change of the stability of 
the fixed point Rq of the deterministic flow f{R) and leads to a signiflcant alteration 
of the ergodic properties. In a vicinity of Rq the behaviour of the stochastic system 
is dominated by the deterministic flow. The transition manifold 

Ac = {a e R'^ : Af = 0} (58) 



does not depend on the system size and the amplitude Ag in contrast to (57). If the 



system lives on the d — 1 dimensional transition manifold (58), that is A/ = 0, the 
scaling of the deterministic flow is not described by (53) but by f{Ro + ~ BfS'^f 



with Uf > nif. The systems behaviour, now depending on Bf,Ag,nf,mg and A^, 
could be classifled in more detail repeating the above procedure. 

4. A^ coupled Stratonovich models 

Now we return to our specific example and consider A^ globally coupled Stratonovich 
models in the strong coupling limit. For drift and noise function we have 

f{R;a) = aR-R^ and g{R;a) = aR. (59) 

The common zero of / and g is Rq = with nif — 2mg = —1, and Af = a, Ag = o. 



Inserting this in (57) we obtain an explicit representation of the transition curve in 
the 2-dimensional parameter space, 

= |(a,a) G : A^- l + 2A^a/a2 = 0} . (60) 

Given the noise strength cr we have 

«c(Ar) = -^fl-l), (61) 



2 V A^. 

which reproduces the results for A^ = 1, for A^ = 2 (see above), and for A^ oo. 



Figure [3] compares results from simulation and the asymptotic result (61) for ac{N) 
and illustrates the finite size scahng adN) — 0^(00) = a'^/(2N) for strong coupling 
> 1. 
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For initial values Xj > (or Xj < 0) \/i the stationary distribution for the center 
of mass (51) lives on 5+ or S-, respectively, and is given by 

( 6{R) if a < a^iN), 

Z = ^((TVAr)i^("-"=(^))r((a-ae(A^))A^/(T') . (63) 

For initial values Xj = Vi we have Ps(-R) = 6{R) for all values of a. 

Similar to the single Stratonovich model, there is a qualitative change in 
the shape of the spatially extended probability density. The maximum of Ps{R) 
undergoes a bifurcation at a™^^ = a'^{l/N — 1/2). Figure lil compares for different 



system sizes the asymptotic result (62) with histograms obtained by simulations for 
large D. For a > a^N) the nth moment of the center of mass can be evaluated as 



a 



f T{{a-a,{N))N/a^ + n/2) 



.N ) r((a - a,(Ar))iV/cT2) 
Keeping N finite, the first moment scales as a — > ac(A^) like 

(i?)^ ~ ± — A(a - a,{N))\ /5 = 1, (65) 




Figure 3: Transition point vs system size N for cP' = 0.1. Simulation results for Z) = 81 obtained 
by maximizing the linear correlation coefficient [22] (gray bullets) and by short time relaxation |23] 
(black bullets) are both in good agreement with the asymptotic result (61 1 for Z? — > oo (solid line). 
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Figure 4: Probability density of the center of mass obtained by simulation for D = 100 (symbols) 
compared with analytic results (lines) for D ^ oo given by Eq. ( [62| ). On the left we show results 
for different system sizes N = 2 (circles), 4 (squares), and 8 (triangles) for a = 1 > a™^^. The 
histograms are obtained from 10^ realizations generated by a stochastic Euler method [26]. On the 
right we show only results for = 8 for Oc < a = —0.42 < a™*'^; here 5 x 10^ realizations were 
generated by a stochastic Runge-Kutta algorithm [l^. — 1. 

since T{z) ~l/zasz— i>0 (25 • Note that also the higher moments (i?") scale hnear 
with a — ac{N). 

Keeping a finite distance to ac{N) we obtain for ^ oo, observing T{z + 
l/2)/T{z) = v^(l - 1/(8^) + . . .) as ;z ^ cx) [25], 

~ ±(a - ae(oo))^ /3 = 1/2, (66) 

which reproduces the result in |2] for D > a^. Higher moments of order n scale with 
(3 = n/2. 

We define the crossover value a*(A^) by (a*(A^) — ac{N))N / a"^ = 1. For 
a <^ a^{N) = —(7^/2 + {3/2)a'^/N we have a linear scaling as for = 1 whereas for 
a 3> a*(A^) a square root behaviour as for ^ oo is observed, cf. Fig. [5] 

Our results are analytically derived for the strong coupling limit in a controllable 
approach. We note that both the critical and the crossover value of the control 
parameter are in accordance with the values proposed on different grounds in [10] for 
weak and intermediate noise, provided the shift due to the Ito interpretation used 
there is taken into account. 
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0.001 0.01 0.1 1 10 100 

a — ac{N) 



Figure 5: Crossover in the scaling behaviour of (i?) as a function of a — adN) as predicted by Eq. 
(64 1 for iV = 1 (dash-dotted line), 2 (dashed line), and 100 (solid line); = 1. The symbols show 
averages over 2 x 10^ realizations generated by a stochastic Runge-Kutta scheme |16] for N = 2 
(circles) and 100 (squares); — 1 and D — 100. The dotted straight lines have slope 1 (left) and 
slope 1/2 (right), respectively. The arrows indicate the crossover points a^,{N) — adN). 



5. Conclusions 



In this paper we have determined the characteristics of a continuous nonequihbrium 
phase transition in a finite array of globally coupled Stratonovich models in the 
limit of strong coupling D ^ oo. In this limit there is a clear separation of the 
time scales governing the evolution of the center of mass coordinate and the relative 
coordinates: The characteristic time of the relative coordinate scales as 1/-D and 
thus becomes short in the strong coupling hmit. The slow center of mass coordinate 
enslaves the fast relative coordinates, its mean value serves as order parameter. This 
allows a controllable and consistent treatment both in the Fokker-Planck and the 
Langevin description which is corroborated by numerical simulations. The reduction 
of a high-dimensional problem to a problem of low dimension is of course inspired 
by generalizations of slaving and adiabatic elimination techniques and the concept of 
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center manifolds to stochastic systems developed in a different context [2211301121], 
cf. also |3I1[32]. 

We have analytically determined both the critical value of the control parameter 
ac{N) and the scaling behaviour of the order parameter and of higher moments. With 
increasing distance from ac{N) a crossover from hnear to square root behaviour is 
found. For N ^ oo the known scaling behaviour is reproduced. The formal results, 
i.e., the computation of the stationary distribution of the center of mass coordinate 
(up to a quadrature) and the determination of the transition manifold are given for 
a general class of systems. 

Our approach may serve as a starting point to calculate next order corrections 
in 1/D. In a multiscale analysis we have to take into account that for finite but large 
D the distribution of the relative coordinates is, though very sharp, of finite width. 

The observation that a solution of the stationary Fokker-Planck equation which 
is not normalizable in a naive way converges to the weak solution if weakly normalized 
is certainly of interest in a broader context. 



Appendix A. Weak normalization 

The FPE for multiplicative noise has two types of stationary solutions: weak 
solutions, i.e. Dirac-distributions living on the zeros of the stochastic flow and 
spatially extended strong solutions which live on a support which is bounded by 
zeros of the stochastic flow or by natural boundaries at infinity. Under certain 
conditions the spatially extended solution may diverge at a zero of the stochastic 
fiow so strongly that it is not normalizable and therefore cannot be considered as 
a probability density. Here we introduce the concept of weak normalization and 
show that in the latter case the weakly normalized solution converges to the Dirac 
distribution living on that zero. 

We assume that the unnormalized solution Ps{x) lives on [xo,b) where Xq is a 
zero of the stochastic fiow and scales for 

Ps(x) ~ const (x - xo)", a < -1 . (A.l) 

The normalization integral diverges at the lower boundary and scales like 

~ 1 
/ dxPs{x) ~ -const A"+^ as A ^ 0. (A.2) 
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Introducing a test function ^{x) which can be expanded near xq as ip{x) = ip{xo) + 
ip'{xo){x — xo) + . . . we have as A ^ 

-b 



r ~ 1 
/ dxPs{x)'^{x) const A"+^ 



a + 1 

X {<^(a;o) + — r<^'(a;o)A + ...}. (A.3) 
a + z 

Now we obtain for the hereby defined weakly normahzed probability density Pf{x) 

!l^/^dxPs{x)^{x 

L 



r ^.dxPsixjipix) 

/ dxPrixMx) = lim ' 7^ = y(xo), (A.4) 

which imphes that P^{x) = 6{x — Xq). 

Appendix B. Langevin approach 

In Sections [2] and [3] we used the Fokker-Planck approach in the center of mass 
and relative coordinates to calculate a^N) for D ^ oo. In this limit the relative 
coordinates — > 0, and it is easy to calculate the reduced stationary probability 
density of the center of mass coordinate Ps(-R). We determined adN) such that 
Ps{R) is a Dirac measure at i? = for a < ac{N) and it is spatially extended for 
a > ac{N). 

The same result can be obtained in the Langevin approach, both in Stratonovich 
and Ito-interpretation as explained in the following for the special case N = 2. The 
generalization to > 2 is straightforward. 

We exploit that for large D the characteristic time scale of the relative coordinate 
r{t) is 1/D so that r{t) becomes very fast. Then the (slow) center of mass coordinate 
R{t) feels only the average of the fast process r(t) and it is justified to replace in 



the Stratonovich-Langevin equation (14) for R the terms associated with r by their 
average, 

dR = {aR - R^ - 3R{r^)) dt+^ (^RodWi{t) + {r o dmit))^ (B.l) 

since for D ^ cxd, r — > we have (r^) = 0. However, the second average 
{r{t)odW2{t)) is not zero as one could naively expect. With the help of the Furutsu- 
Novikov theorem |33l EH we obtain 



(r(t)o6(t))= / dsums)){P^^ 
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1 / 6r{t) 



2 \SUs) 



1 a 



{Rit)) . 



(B.2) 



Note that the averages here are with respect to reahzations of ^2- 

We now observe that the resulting equation for R does not depend on ^2, 
therefore {R) = R, and obtain 



dR 



a + 



R-R-" 



dt+^RodWi{t) 
v2 



from which the threshold ac(2) = — follows. 

The system ( I4p5 ) in Stratonovich sense can be written in a compact form as 
dp = i{p)dt + J2j=i 2 (p) ° '^^i(^)? where p = {R, r)^. The equivalent Ito system 
is dp = (f + 1/2 J2j dpg'^^^g^^^)dt + J2j S^^^ dWj{t), where the drift term is modified 
by the Ito shift; dpg'^^^ is the shorthand of a Jacobian, cf. e.g. For our system 
we have g*^^^ = a/\/2 p and g*^^^ = (t/a/2 {r,R)'^. The Ito shift amounts to cr^/2 p 
so that the equivalent Ito version of (14) reads 

,2 



dR 



a + 



a 



R-R'-mr' 



a 



dt + — [RdWi{t)+rdW2{t)y (B.4) 

Again, R feels only the average of the terms associated with the fast process r, we 
have (r^) = but now also the second average vanishes since in the Ito calculus 
{r{t)dW2{t)) = {r{t)){dW2it)) = which results in 



dR 



a + 



a 



R-R-" 



dt+ ^RdWiit). 
v2 



This is indeed the Ito eqivalent to Eq. (|B.3|) which can be seen observing that in the 

- l/N) as 



single variable case the Ito shift is simply 1/2 g' g = a"^ /A R. 

For arbitrary the same procedure leads to ac{N) = — ((t^/2) 
obtained above. 
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